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(57) Abstract 



Complete helical cone-beam scanning and non-redundant 
data acquisition are obtained for three-dimensional tomographic 
imaging of infinitely long objects. The minimum sized 
two-dimensional detector window is bounded by two consecutive 
turns of the helix. The ray source exposes all object points 
during a rotation of exactly 180 degrees when seen from the 
points themselves. Only one-dimensional filtering is employed 
in the reconstruction. Rebinning to parallel beams, as seen along 
the axis of rotation, allows for especially simple procedures 
without any need for pre-weighting or magnification factors. As 
a special case, the invention is applicable to helical fan-beam 
scanning with one-dimensional detector arrays. 
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1 

Technique and arrangement for tomographic imaging 
Background 

5 This invention was first described in [Dan97a]. A two-dimensional detector 11 and a 
point-shaped ray source (e g X-ray) S are assumed to move synchronously around the object 
in a helical trajectory as shown in Figure 1. In a medical tomograph the helical source 
movement is achieved by translating the patient through the rotating source-detector gantry 
with constant speed. Two-dimensional projections are acquired (detected) at arbitrarily short 

10 intervals along the trajectory. The detector consists of a large number of sensors (detector 
elements) which are evenly spaced and placed in a plane or, as in Figure 1, on the surface of 
the helix cylinder 12. Although the rotation axis 14 is normally horizontal in medical 
tomographs, rather than vertical as in Figure 1, we will adopt the following convention. In 
the sequel, vertical means a direction parallel to the rotation axis 14 (the z-axis) in Figure 1, 

15 while horizontal means a direction parallel to the xy-plane 15. 

A projection consists primarily of intensity measures for the incoming rays to a detector 
element The logarithms of these primary data represent the sum of the attenuation along the 
rays, i.e. line integrals over the three-dimensional attenuation function / we want to retrieve. 

20 But to be able to reconstruct / from its projections in a correct way, all points in the object 
have to be fully exposed and the projection data utilized in a balanced way. Thus, if back- 
projection is used for the reconstruction, projections from all projection angles must be 
available and brought in with the correct weight to obtain what is called exact reconstruction. 
Also, the projection data have to be filtered correctly to compensate for the inherent low-pass 

25 filtering in the projection-back-projection procedure. 

In the literature several mathematically exact methods have been proposed for reconstruction 
from cone-beam projections. In most cases these methods demand that the object is of finite 
extension, i.e. restricted in size, so that its total projection never falls outside the available 

30 detector. Unfortunately, this requirement is not realistic in most cases of computer 
tomography, e.g. when one is to reconstruct a full body, or long objects in general. 
Traditionally, ID-detector arrays are made large (wide) enough to cover the object across its 
maximum width. However, for several reasons, it is out of question to extend these ID- 
detectors to 2D-detectors, which cover the patient from head-to-toe. Instead, in the 

35 foreseeable future, available 2D-detectors will be used to cover and record projections of a 
section of a long object. 
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Today, three-dimensional volume data are reconstructed slice-by-slice. The patient is 
translated slowly (typically 2 mm/sec) while the X-ray source and a one-dimensional detector 
array are synchronously and continuously rotated at speeds of around 1 r/sec. Relative to a 

5 patient which is not moving, the source and detector are then performing a helical movement 
with very low pitch, say, 2 mm. The reconstruction employs a modified versions of 
traditional 2D-reconstruction methods for circular scanning of a single slice. However, with 
the given numbers, it takes approximately 100 sec to fetch data for a 200 mm long section of 
the body. During this time, due to breathing and other body functions, the body is not fully at 

10 rest which blurs the reconstructed object. A second drawback is that the anode of the X-ray 
tube is subjected to severe strain and extreme temperatures during longer exposure times. 

In a ID-detector system the major part of the generated photons are collimated away without 
being utilized, while a 2D-detector system is able to utilize a substantial part of these 

15 otherwise wasted photons. Hence, by using a 2D-detector with, say, n parallel ID-detectors 
in the above example, the velocity can be increased to 2n mm/sec and the scanning time 
reduced to 100/n sec. Alternatively, speed can be traded for strain on the X-ray source so 
that, for instance, if the photon flow is halved, the velocity is more moderately increased to n 
mm/sec and the scanning time reduced to 200/n sec. However, in any case it is no longer 

20 possible to perform the reconstruction using conventional 2D-methods since the projection 
rays are no longer, not even approximately, in the same plane during one turn of the source 
trajectory. 

25 



30 
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Prior Art 

Circular source trajectory 

5 A well-known method for inexact reconstruction from cone-beam projections taken along a 
circular path was proposed in [Feld84]. The 2D-detector is placed on a planar surface and 
extended horizontally to cover the width of the object. The width of the object and its distance 
to the source defines the maximum fan-angle y max of the source-detector system. In the 
vertical direction the planar detector is limited by two horizontal lines. Along the vertical axis 
10 where these lines are closest to the source we find the maximum cone-angle. The image 
reconstruction consists of the following steps taken for each detector recording. All 
corrections of geometrical and radiometric nature, including the ever necessary logarithm 
computation have been left out here for the sake of brevity. 

15 1. Pre-weighting of the recorded detector data with a factor that is proportional to the cosine ' 
oftlie angle between the central ray and the ray that originated the detected value. 

2. Filtering with traditional ramp-filtering techniques along each horizontal detector row. 

3. Back-projection along the original ray in which process the filtered detector value is 
multiplied with a so called magnification factor which depends on the distance between tlte 

20 ray source and the object point to receive a contribution from the ray. 

This method gives perfect results for image slices in, or close to the mid-section of the object. 
For slices which have been subjected to more oblique rays at higher cone angles, the image 
quality deteriorates. 

25 

Helical source trajectory. Non-exact methods 

Extensions of [Feld84] to helical source paths were first proposed by [Wang93]. Here, the 
planar 2D-detector is given a vertical extension large enough to ascertain that every point is 

30 exposed to the source at least once for every projection angle during a full 360 degree source 
rotation. The effect of this requirement is that for any given projection angle an object point 
will be exposed by the source from various numbers of source positions; at least one 
but often many more, depending on the given fan-angle, cone-angle and detector size. This 
has to be taken into account during the back-projection. Hence, [Feld84] is employed in 

35 [Wang93] but augmented with the following rule. 
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3a. During the back-projection, Jot a certain projection angle, among all possible source 
positions which illuminate an object point, contributions are accepted only from the 
position which is closest to the actual point in tlie z-direction. 

5 A way to achieve a more efficient and balanced exposure of the object points was proposed in 
[Scha96]. Here, the detector is located (wrapped) onto the surface of the source cylinder 
41, which in Figure 4 is seen to be centered in S. The radius of this cylinder equals the 
source-detector distance, which is different from the radius R of the helix cylinder 12. The 
helix cylinder is coaxial with the object cylinder 13 in Figure 1 which is defined by the 

10 maximum object width r. In [Scha96] the detector is limited in the vertical direction by two 
horizontal circles (cross-sections) of the source cylinder 41, However, it is not quite clear 
what the minimum or optimal height is to be recommended for the detector. In the horizontal 
direction the detector is limited by two vertical lines, set to let the detector cover the object 
cylinder. In the following we may use Figure 1 to clarify some prior art such as this, even if 

15 the detector arrangement in our invention is different from the detector in [Scha96] and 
[Scha97]. 

The main novelty in [Scha96] is the introduction of complementary projections. These are 
projection data captured at the source cylinder 41, but sorted and resampled (rebinned) with 

20 respect to where the rays from the source are reaching the helix cylinder 1 2 . Assume for the 
moment that the source is moving along the helix 16 while the object and the helix cylinder is 
fixed. All rays from various source positions on the same turn, reaching a fixed point on the 
helix cylinder, is said to be a complementary fan-beam and the projection data for this fan- 
beam is created during the rebinning. The set of such not quite horizontal and not quite planar 

25 fan-beams, with rays fanning out from points on a vertical line on the helix cylinder, 
constitute a complementary projection. These are employed in the following reconstruction 
procedure. 

1. Pre-weighting of the recorded detector data with a factor that is proportional to the cosine 
30 oftlie angle between tlie central ray in tlie projection and the ray that origuuited ttie 

detected value . 

2. Re-binning to complementary projections. 

3. Filtering of the original as well as the complementary projections with traditional ramp 
filtering techniques along each liorizontal detector row. Because of the non-planar detector 

35 this filter is slightly different from tlie filter employed in [Feld84J. 

4. For each projection angle, filtered projection data are back-projected along the rays. The 
values are multiplied with magnification factors which depend on the distance between tlie 
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point and ray source. All such values, from original as well as complementary corn-beam 
source positions, are averaged into one single contribution that is accumulated to the object 
point. 

5 We notice that the detector arrangement in [Scha96] does not secure a perfectly balanced 
exposure of the object points. During the back-projection event, for each rotation angle, there 
is a similar situation as in [Wang93] where the object points are exposed from one or several 
source positions. The difference is that in [Scha96] all these projection data from both 
original and complementary projections are utilized and averaged together during the 
10 backprojection. 

[Scha97] proposes another reconstruction technique that is claimed to be more 
computationally efficient. The detector system is identical to the one in [Scha96] with two 
horizontal truncating circles on the surface of the source cylinder 41. The reconstruction 
15 consists the following steps. 

1. Rebinning to oblique parallel projections. 

2. Pre-weighting of the recorded detector data with a factor that is proportional to the cosine 
of tlie angle between the ray tliat originated tlie detected value and tlie central ray. 

20 3. Reconstruction of one horizontal slice from generalized projections. The latter can be 
seen as the result of imaginary projection rays running within the horizontal slice. 

3.1. Computation of Fourier domain contributions to this slice for one generalized 
projection at every projection angle. 

3.1.1. Computation of Fourier domain contributions for each "detector" position in one 
25 generalized projection. 

3 .1 J J . Multiplication of projection data (from all source positiom tfiat send oblique rays 
through tlie slice in this position) with a pre-computed set of weights, which are 
Fourier series components derived from an adopted interpolation function. 

3.1.12. Summation of tlie contributions for each Fourier component to obtain a single set, a 
30 truncated Fourier transform for each ray in this detector position of this generalized 

parallel projection. 

3.1.2. Computation of the Fourier transform (FFT) along the projection for all these 
truncated Fourier components to obtain a kind of truncated 2D Fourier transform 
contribution for each generalized parallel projection. 

35 3.1.3. Multiplication of the Fourier transform of this generalized projection with a ramp 
filter. 

3.2. Merging of filtered data from all projection angles in the 2D Fourier space of tlie 
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liorizontal slice and resampling with a space-invariant interpolation filter. 
3.3. Application of an inverse 2D Fourier transform (FFT). 

3 A. Compensating for imperfect interpolation in tlie Fourier domain by post-weighting the 
result with tlie inverse interpolation function, in accordance with the well known 
5 gridding technique. 

The first rebinning step is best understood if the source, the detector, and the object 17 is 
pictured as seen from above. From there, the cone-beams will be seen as fan-beams. The 
rebinning in step 1 above is equivalent to sorting projection data into sets where data from 
10 this point of view are produced not by fan-beams but from parallel beams. The term oblique 
parallel projections stems from the fact that the rays are parallel when seen from above, but 
in general non-parallel and oblique to the horizontal plane. To understand the following steps 
it is now recommendable to imagine a planar, virtual detector 122 as in Figure 12 placed 
vertically on the rotation axis in the middle of the object. There are several source positions 
15 which produce the rays for this projection. Since the real detector on 41 is truncated 
horizontally and the source positions are located along the helix 16, the effective area of this 
virtual detector does not have a left-right symmetric shape. The upper and lower boundaries 
131 and 132 are curved and tilted as shown in Figure 13. This is an important difference to 
the perfectly rectangular shape of the corresponding virtual detector 72 in Figure 10 for the 
20 present invention. The net effect is that in [Scha97] a varying number of source positions 
generate fan-beams which penetrate a given slice under various oblique angles. All of these 
contribute to the result in the above ingenious but rather complicated computation steps 
3.1.1.1 and 3.1. 1.2. In the present invention we will find one and only one such fan-beam. 

25 Helical source trajectory. Exact methods 

An exact method for reconstruction of a limited sector, a Region Of Interest (ROI) of a 
long object was proposed in [Tam95] and [Eber95]. The helical scanning covers the full 
vertical extension of the ROI but has to be complemented with two circular scans at top and 

30 bottom, respectively. The detector is placed on a planar surface, just as in [Feld83] and 
[Wang93] but the detector window is limited to the area between two consecutive turns of 
the source trajectory 16 as in Figure 1. The upper and lower truncating lines on the detector 
plane are therefore neither horizontal, straight, nor left-right symmetric. The arguments for 
the specific extension of the detector stem from a well-known completeness condition for 

35 Radon planes which carries over to the following reconstruction technique. In essential 
aspects this method is an outgrowth from [Gra87]. 
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1. From each 2D-projection, partial contributions to the derivative of Radon transform values 
are computed by tneans of line integration along a multitude of Urns in the planar detector. 
This requires that we select a specific object point to be the origin of a 3D coordinate 
system. 

5 2. When the scanning is complete, that is when the helical trajectory has covered the intended 
target region oftlie object (ROI), all tliese partial contributions are sorted and coplanar 
partial contributions are summed. 
3. The result is resampled into a regular grid in the Radon transform space of the R01 of the 
object function. 
10 4. Filtering with a derivative filter. 

5. 3D back-projection which takes place as two consecutive 2D back-projection steps. 

Generally speaking, this reconstruction is more complicated and costly that the previous 
ones. Also, the rhythm of the reconstruction procedure is affected by the chosen size of the 
15 ROI. It does not feature an even flow of identical procedures repeatedly taking place for every 
new projection regardless of the length of the object. The two extra circular scans are highly 
unwanted since they break the smooth and continuous translation-rotation motion of the 
helical part. However, the method is optimal in one respect. For a given pitch of the helix it 
utilizes a minimum sized detector. 

20 



25 



30 
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The Invention 

The present invention utilizes the same optimal, minimum cost two-dimensional detector 
geometry as proposed in [Tam95], characterized by an exposure window which is limited 

5 vertically by the two nearest turns of the helical source trajectory. However, both the 
motivation for and the exploitation of this detector window differs greatly from the ones 
given in [Tam95] and [Eber95]. To explain the specific virtue of this exposure window, we 
refer again to Figure 1, which shows a perspective view of a source S, a detector 11 
wrapped around the helix cylinder 12, and inside this an object cylinder 13. In the sequel, 

10 unless stated otherwise, we assume that the object cylinder is rotating counter-clock-wise as 
shown around thez-axis and translated upwards in a right-handed helix, while the source S 
and the detector 11 are fixed in the space (x, y, z). 

Figure 2 shows the arrangement as seen from above, while Figure 3 shows the detector 
15 window unwrapped and rolled out on a plane. Note that Figures 1 and 3 are consistent only if 
the rays in Figure 3 are understood to be coming from the source towards the viewer. Figure 
4 shows the detector placed on the source cylinder 41 centered in S and having a radius 
which is twice as large as the helix cylinder 12. The figures 1, 2, 3, and 4 presented here are 
in principle, although not in detail, identical to figures 1, 2, 3, and 4 in [Dan97a]. 

20 

As mentioned, the 2D-detector 11 in Figure 1 is wrapped onto the helix cylinder 12. 
Unwrapped and rolled out on the plane of the sheet, the same detector surface 1 1 in Figure 3 
is seen to be bounded by four straight lines, two vertical ones 31 and 32, and two slanted 
ones 33 and 34. Within this area the object 17 is projected, i.e., rays from the cone-beam 
25 source reaches active detector elements. Horizontally, this area has to be extended to cover 
the object cylinder 13, which translates to a certain width, or fan angle , as seen from 

the source. As an example we have assumed that this object cylinder has a radius r = -j= 

where R is the radius of the helix cylinder 12. This means that horizontally on 12, the 
detector covers a rotation angle of 180 degrees out of 360, and that seen from the source the 
30 detector 11 covers a fan angle from - 45 to +45 degrees. In principle the detector may be 
extended to a full turn which then has a fan-angle from -90 to +90 degrees and would allow 
for an object cylinder that extends all the way to the helix. The slanted lines 33 and 34 are 
intersecting the cylinder surface 12 at the slope 

tane = — « (1) 
coR 2jtR 
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where v is the vertical translation velocity, <d is the angular velocity for the rotation, and h 
is the pitch of the helix. 

At the core of the invention is the following property of the detector-exposure window. 

5 Every point in a cylindrical long object, with a radius that fits inside the 
boundaries of the detector window, will be exposed (projected) during a 
rotation angle which is exactly 180 degrees, seen from the actual point in the 
object. A more formal proof of this property was given in [Dan97b] but is left out here for 
the sake of brevity. A conjecture of this new sufficiency condition is that as soon as one 

10 point or a set of points (i.e. a part of the long object) has been fully exposed 
in the above sense, the reconstruction of this part can take place. This is in 
contradiction to the situation in |Tam95] and [Eber95]] where the whole ROI has to be 
exposed to make the Radon space complete before the actual reconstruction is commenced. 

15 An example of this 180 degree exposure is the line 18 in Figure 1. It contains the three object- 
points Q x - Q - Q 2 and it is shown in two positions where the exposure starts and ends, 
respectively. Note that the end points of this line is sliding and touching the outer cylinder so 
that during the rotation, both ends will coincide with the source S. Any such line will be 

called a Jt-Iine. This line is also shown in Figure 2 in the same two positions. 

20 

Assume as before that the object is moving upwards and rotating counter-clockwise when 
seen from above. In the detector window of Figure 3 the line Q x -Q-Q 2 crosses the lower 
boundary 34 as a single point at Q Ul . After a rotation with the angle n + 2y around the axis 
14, this line will be seen as a single point again from the source leaving the detector at Q OM 
25 on the upper boundary 33. Clearly, between entrance and exit the source has rotated exactly 
180 degrees as seen from any point on this line. Since we have chosen this line quite 

arbitrarily, the same thing is true for all points in the object which belong to fully exposed %- 

lines. In Figure 2 and Figure 3, but not in Figure 1, we have inserted another Jt-line 

P x -P-P 2 . In the fixed source-detector system of Figure 2 this line P enters and exits in 
30 positions which are exactly the reverse of the corresponding positions for the line Q. The line 
P is therefore closer to the source than line Q during its exposure, which takes place during a 
rotation angle of n - 2y . The points on line P travels over the detector surface along different 
and shorter curves as shown in Figure 3, but seen from any of these points, the source 
rotates around them exactly 180 degrees. 

35 
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It was shown in [Dan97b] that every object point belongs to one and only one irc-line. 

Therefore, the detector system in Figure 1 gives us a complete and perfectly balanced data 
capture for every point and hence also for the whole object. Furthermore, from the conjecture 
above follows that it should be possible to reconstruct the object at the same pace as an 

5 incremental part (each new set of rc-lines) of the long object is fully exposed. 

The physical implementation and placement of the detector can of course be made in various 
ways as indicated in Figure 4. For instance, it may be placed on the helix cylinder 12 itself, 
on the source cylinder 41 or on a plane 42. In any case, the detected and utilized data must 
10 be restricted to the window defined by Figure 3. 

In our invention, using the same detector data, the elaborate reconstruction in [Tam95] and 
[Eber95] will be replaced by a much simpler procedure. To describe this procedure, we do 
not have to limit the ongoing scanning and reconstruction to a predetermined ROI, nor do we 
15 have to specify a 3D origin for the process. Instead, scanning and reconstruction is like a 
constantly ongoing flow, in principle without beginning or end, where each new projection is 
absorbed and incorporated seamless to the previous result. For this purpose, the following 
general reconstruction procedure for every new projection is found in [Dan97a]. 

20 l.Rebinning 

2. Pre-weighting (depending on rebinning and detector type) 

3. One-dimensional filtering with a ramp-filter across the detector {where the filter design is 
dependent on rebinning and detector type) 

4. Back-projection along incoming ray direction with magnification factors, depending on 
25 type of rebinning as well as on detector type: plane, cylindrical etc. 

A special case of this procedure is rebinning to parallel projections which we will describe in 
more detail. Figure 5 shows a straight side-view of Figure 1 with six rays 51, 52, 53, 54, 
55, and 56 coming from the source S positioned at the x-axis. Figure 6 shows a view from 

30 above where the object is fixed and the source and detector is rotating. With the source in the 
position S a we observe three fan-beams 61, 62, and 63 (seen as rays in this view), which 
comprises the six rays in Figure 5 and which produce the three projection sets 
/(a, y,), *(a,.0), /(a,-y,). The two outer rays are parallel to two other rays, 64 and 65, 
coming from two other source positions which produce the projections r(a+yj, 0) and 

35 t(a -y p 0), respectively. Clearly, we may resample our projection data so that data from 
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such parallel fan-beams (seen as rays) are brought together. This can be done with either of 
. the following two equivalent assignments. 

[/*a+y,y) <= *(a,y)] - [rffty) *(/3-y,y)] (2) 

5 

As shown in Figure 7 we are then free to see the data set p(/3, y) as generated by a parallel 

beam in the /J-direction. Without loss of generality, this direction is horizontal in Figure 7. 
Perpendicular to these rays we place a virtual detector 72 on a vertical plane. 

10 The detector window 71 for the parallel projection in Figure 7 is unwrapped and rolled out 
into the sheet of Figure 8. Note that the complete detector positions for the parallel projection 
are put together from vertical lines 83, 84, and 85 each one stemming from different cone- 
beam detector positions. The resulting parallel beam detector area has the same slant as the 

cone-beam detector but is shortened with a factor of two in the y-direction. The uppermost 

15 and lowermost part of the detector 81 and 82 in Figure 8 outlines another detector window 
included here for comparison only. To the best of our understanding, this window 
corresponds to the minimum size detector in [Scha96] and [Scha97] when mapped onto the 
helix cylinder 12. For the given pitch = h and the given maximum fan angle y^ the height 
of this detector window is 



2Q 2v( J r+y niax ) m h x + Ym * m h 

(oR cos y msix cosy,™ ji 



\ 4. ^nrax \ 
71 



cosy r 



(3) 



This formula indicates that the detector redundancy in [Scha96] and [Scha97] grows rather 
quickly for increasing fan-angles. 

25 The rays in the parallel projection emanate from a set of sources with vertical fan- 
bcams, located on a specific section of the helix. Rolled out in the plane of the sheet this part 
73 of the source helix is superimposed on the detector 71 in Figure 8. It takes the form of a 
line with the same slant as the detector but with opposite sign. Because of this fact, in the 
present invention, the virtual detector 72 in the vertical mid- plane is bounded by a perfect 

30 rectangle with a width that equals the object cylinder diameter and a height which is exactly 
half the pitch = /z/2. This is illustrated in Figure 9, where an upward tilt of the source path 7 3 
is exactly compensated for by a downward tilt of the detector. Furthermore, since the distance 
from the virtual detector 72 to the source is everywhere identical to the distance to the real 
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detector, the real detector height h is always demagnified to exactly hJ2 at the virtual detector. 
Figure 10 illustrates the second part of the rebinning-resampling procedure, namely from 

equidistant grid points in y R to equidistant grid points in y 1 = R siny . 

5 The aforementioned property of the virtual detector area being a perfect rectangle is further 
illustrated in Figure 11, which shows three orthogonal views A, B, and C of the parallel 
. projection system. Seven source positions are indicated. In A, B we can see the projection as 
a line d-e from the source positions S. Clearly, in view B we see that all the three points S , 
d, and e are on the helix. Furthermore, the plane of the virtual detector intersects the helix in 

10 two points which are exactly halfway between S and d at the upper ray 11 1 and halfway 
between S and e at the lower ray 112. Therefore the height of the vertical detector is ft/2 
with its midpoint on the x-axis for any S. This proofs that the virtual detector is a rectangle 
with horizontal boundaries. 

15 Thus, using the insight that there is a special detector window which delivers sufficient and 
non-redundant data, we capture cone-beam projection data on this detector and rebin them 
into parallel projection data to create an advantageous situation for the actual reconstruction. 
The complete procedure consists of the following three steps. 

20 1. Rebinning to parallel projections as described by tlie Figures 6, 7, 8, 9, 10 t and 11. 

2. Filtering with a conventional ramp- filter along liorizontal rows in the virtual detector plane. 

3. Back-projection in the direction of the original rays using a constant magnification factor. 

In the present invention, after parallel rebinning, the one-dimensional filtering takes place 
25 along horizontal rows in the virtual detector 72 of Figure 10. In [Scha96] and [Scha97] the 
filtering takes place along horizontal rows of a detector placed on the source cylinder arc 4 1 
in Figure 4. Figure 13 shows this detector mapped onto the virtual detector plane 121. The 
horizontal rows in the real detector are mapped onto curves in 121 which are neither 
horizontal nor straight. Clearly, filtering along such curves in the virtual detector plane rather 
30 than along horizontal rows will produce a rather different reconstruction result.. 

Even so, Step 3 in the above procedure may very well be replaced by a simplified version of 
the corresponding step 3 in [Scha97]: Reconstruction of one horizontal slice from generalized 
projections. The simplification is due to the perfectly balanced data capture in the present 
35 invention. We know a priori that there is one and only one source position that contributes to 
each detector position in the generalized projections as shown in Figure 9. Hence, there is no 
need to keep track of multiple exposure contributions, since there are neither missing nor 
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redundant data in any projection. The situation is different in [Scha97] which is illustrated in 
Figure 12 showing a vertical section of the parallel scanning system. The real detector 125 is 
much higher than in Figure 9 so that the virtual detectors 121, 122, and 123 for 
neighboring half turns overlap vertically. Therefore, in a vertical plane (such as the plane of 
5 the sheet) a horizontal slice of the object is partially illuminated not from one but from three 
source positions on the trajectory. This irregularly distributed redundancy in exposure is also 
reflected in Figure 13 which shows the virtual detector window in [Scha97] for the minimum 
sized detector. The upper and lower boundaries 131 and 132, respectively, are the same as 
81, 82 in Figure 8, although mapped onto the virtual planar detector. 

10 

In the most likely physical embodiment of the 2D-detector arrangement proposed in this 
invention, the detector elements are placed onto the source cylinder 41. See Figure 4. For 
moderate cone angles the detector elements are then facing the incoming rays rather straight 
on. For detectors made to cover high cone angles it might be more appropriate to mount the 
15 detector elements on the inside of a sphere centered in S. This would guarantee or at least 
make it more easy to secure that all detectors are facing the incoming rays correctly. 

Figure 14 shows again the detector window 11 on the helix cylinder rolled out on the plane 
of the sheet. However, this time it is overlaid with the same the detector window mapped 
20 onto the source cylinder arc 41. In its rolled out version, this detector appears in Figure 14 
outlined as 141. Considering the geometry of Figure 4, it might be more optimal to place the 
detector on the source cylinder arc 43 having the smallest possible radius close up to the 
object cylinder 13. However, since the geometry of such a detector would conform exactly 
with 141, we may discuss the geometry of 141 without loss of generality. 

25 

The detector 141 coincides with 11 in the middle but varies with y so that the top-most . and 
bottom-most point of the detector are found at 



z top= — and z bottom = — , (4) 

n J? rr\Q v ml? mc v 



** coR cosy bottom ujR cosy 

30 

respectively. The height H is then varying as 



v Tji 

cuR cosy cosy 



= Z lap + Zbouom " 37 7T~ - 7Z - ( 5 ) 
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where h is the pitch as before. Thus, data which are captured on the source cylinder have to 
be resampled from the unevenly sloping detector area in Figure 14 to the grid of the detector 
in Figure 3 (also shown in Figure 14), defined by vertical lines and evenly sloping lines with 
rhombus shaped detector elements. When projection data are resampled once more to the 
5 planar virtual detector in Figure 10, the final grid pattern will be perfectly rectangular. 

An important special case for the present invention is when the detector 141 (and the pitch) 
of Figure 14 is reduced in height to a single row 150 of detector elements 151, which is 
shown in Figure 15. We note that also in this special case will the height of the detector 

10 element increase with increasing fan angle as predicted by the above formula (5). Normally, 
the detector array in Figure 15 would no longer be considered as a two-dimensional detector 
but a one-dimensional array detector. One-dimensional array detectors are used in existing 
helical fan-beam tomographs for which the state-of-the-art is represented by [King93], The 
detector is normally placed on the surface of a source cylinder 41, although not designed as 

15 the one in Figure 15. Instead, the detector elements are of constant height and they are not 
placed in a slanted fashion but horizontally straight on the source cylinder surface. 

As a consequence, to secure sufficient data, either the height of the detector elements have to 
be increased, as in formula (3) which decreases the resolution in the z-direction, or the pitch 

20 of the helix has to be decreased with the same factor, which reduces the scanning efficiency 
and increases the dose compared to the present invention. The scanning will also acquire 
much redundant data so that the accompanying reconstruction procedure has to employ 
elaborate weighting factors to compensate for multiple exposure. In essence, [King93] is a 
predecessor and the fan-beam counterpart to [Scha96]. Therefore, in helical fan-beam 

25 scanning, the present invention should have the same advantages compared to [King93] as 
we have found it to have when compared to [Scha96] and [Scha97]. Using the present 
invention with a detector designed and arranged accordingly, for instance as in Figure 15, the 
data capture will be complete and free of redundancy and the reconstruction procedure can be 
simplified to contain the three steps rebinning, one-dimensional ramp filtering, and 

30 backprojection with constant magnification factor. 



35 
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CLAIMS 

1. Arrangement and method for three-dimensional tomographic imaging of long objects, 
where the object to be investigated is subjected to simultaneous translation and rotation 
5 relative to the cone-beam ray source and the detector, characterized by a two-dimensional 
detector limited to a window opposite to the source, where this window is constrained to 
correspond to a surface area shaped as a parallelogram which is attached to a cylinder 
centered at the axis of rotation and passing through the source, and where the height of the 
window measured along the length of the cylinder equals the pitch of the helix and the width 
10 of the window is made to cover the maximum width of the object to be reconstructed, and 
where the upper and lower boundary of the window coincides with two consecutive turns of 
the source path relative to a fixed object, and where the reconstruction procedure is 
characterized by steps which include pre- weighting with a factor that depends on the angle of 
the ray that gave rise to the specific detector value, followed by filtering with a filter of type 
15 ramp-filter horizontally or near horizontally across the detector, followed by back-projection 
with a magnification factor along the direction of the rays that gave rise to the original detector 
values, and where the actual pre-weighting, the actual filter design, and the actual 
magnification factor depend on the physical embodiment given to the detector and which 
rebinning that has been employed to the detected data. 

20 

2. Arrangement and method for three-dimensional tomographic imaging of long objects, 
where the object to be investigated is subjected to simultaneous translation and rotation 
relative to the cone-beam ray source and the detector, characterized by a two-dimensional 
detector limited to a window opposite to the source, where this window is constrained to 
25 correspond to a surface area shaped as a parallelogram which is attached to a cylinder 
centered at the axis of rotation and passing through the source, and where the height of the 
window measured along the length of the cylinder equals the pitch of the helix and the width 
of the window is made to cover the maximum width of the object to be reconstructed, and 
where the upper and lower boundary of the window coincides with two consecutive turns of 
30 the source path relative to a fixed object, and where the reconstruction procedure is 
characterized by steps which for each incoming cone-beam projection includes rebinning to 
parallel projections as seen along the rotation axis and locating and resampling the rebinned 
data to a virtual detector plane on the rotation axis, followed by ramp-filtering along all 
horizontal rows in this detector plane, followed by back-projection in the direction of the 
35 original rays using a constant magnification factor. 

3, Arrangement and method for three-dimensional tomographic imaging of long objects, 
where the object to be investigated is subjected to simultaneous translation and rotation 
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relative to the cone-beam ray source and the detector, characterized by a two-dimensional 
detector limited to a window opposite to the source, where this window is constrained to 
correspond to a surface area shaped as a parallelogram which is attached to a cylinder 
centered at the axis of rotation and passing through the source, and where the height of the 

5 window measured along the length of the cylinder equals the pitch of the helix and the width 
of the window is made to cover the maximum width of the object to be reconstructed, and 
where the upper and lower boundary of the window coincides with two consecutive turns of 
the source path relative to a fixed object, and where the reconstruction procedure is 
characterized by steps which for each incoming cone-beam projection includes rebinning to 

10 parallel projections as seen along the rotation axis and locating and resampling the rebinned 
data to a virtual detector plane in the rotation axis, followed by reconstruction of one 
horizontal slice at a time using generalized projections and Fourier transform technique. 

4. Arrangement and method for three-dimensional tomographic imaging of long objects as in 
15 Claim 1, 2, or 3, characterized by an effective detector area that has an extension which 

correspond to the area exposed through the aforementioned window but is physically shaped 
in a different way and placed in a different position. 

5. Arrangement and method for three-dimensional tomographic imaging of long objects as in 
20 Claim 4, characterized by a detector which is placed on a vertical plane. 

6. Arrangement and method for three-dimensional tomographic imaging of long objects as in 
Claim 4, characterized by a detector which is placed on the surface of a cylinder with a 
vertical axis on the source S. 

25 

7. Arrangement and method for three-dimensional tomographic imaging of long objects as in 
Claim 4, characterized by a detector which is placed on the surface of a sphere centered on the 
source S. 

30 8. Arrangement and method for three-dimensional tomographic imaging of long objects as in 
Claim 4, characterized by a detector which is placed on the surface of a cylinder, the axis of 
which is the axis of rotation. 

9. Arrangement and method for three-dimensional tomographic imaging of long objects as in 
35 Claim 4, 5, 6, 7, or 8, characterized by a detector consisting of one single row of detector 
elements. 
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Figure 2 
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Figure 4 
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Figure 6 
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Figure 7 
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Figure 9 
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Figure 11 
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